O 
O 
(N 



< 



> 
in 



J2 

Oh 



O 

CO 



To Appear in Proceedings of the Spitzer Science Center 2005 Conference: 
Infrared Diagnostics of Galaxy Evolution 
ASP Conference Series, Vol. tbd., 2006 
Ed. R. Chary 

Spitzer Observations of Red Galaxies at High Redshifts 



Casey Papovich, for the GOODS and MIPS GTO teams 
Steward Observatory, 933 N. Cherry Ave., Tucson, AZ 85721 



Abstract. I discuss constraints on star-formation and AGN in massive, red 
galaxies at z~l-3 using Spitzer observations at 3-24 /mi. In particular I focus 
on a sample of distant red galaxies (DRGs) with J— K s >2.3 in the southern 
Great Observatories Origins Deep Survey (GOODS-S) field. The DRGs have 
typical stellar masses A^IO 11 Mq. Interestingly, the majority (j>50%) of these 
objects have 24 /im flux densities ;>50 /iJy. At these redshifts massive galaxies 
undergo intense (and possibly frequent) IR-active phases, which is in constrast 
to lower-redshift massive galaxies. If the 24 /im emission in these z~l-3 galaxies 
is attributed to star formation, then it implies star-formation rates (SFRs) in 
excess of ~100 M.® yr _1 . These galaxies have specific SFRs equal to or exceeding 
the global average value at that epoch. Thus, this is an active period in their 
assembly. Based on their X-ray luminosities and near-IR colors, as many as 
25% of the massive galaxies at z;> 1.5 host AGN, suggesting that the growth of 
^h ■ supermassive black holes coincides with massive-galaxy assembly. 
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1. Introduction 



Most (~50%) of the stellar mass in galaxies today formed during the short time 
between z~3 and 1 (e.g., Dickinson et al. 2003, Rudnick et al. 2003). Much 
of this stellar mass density resides in massive galaxies, which appear at epochs 
prior to z~l-2 (see, e.g., McCarthy 2004, Renzini 2006). It is still unclear when 
and where the stars in these galaxies formed. The fashionable scenario is that 
galaxies "downsize" , with massive galaxies forming most of their stars in their 
current configuration at early cosmological times, with less-massive galaxies 
continuing to form stars to the present (e.g., Bauer et al. 2005, Juneau et al. 
2005). But it is unclear when and where the stellar mass in galaxies forms. For 
example, it may be that stars form predominantly in low-mass galaxies at high 
redshifts, which then merge over time to form large, massive galaxies at more 
recent times (e.g., Kauffmann &: Chariot 1998, Cimatti et al. 2002). 

At z<^l most massive galaxies exist on a fairly prominent red sequence 
(e.g., Blanton et al. 2003, Bell et al. 2004), are largely devoid of star formation, 
evolve passively, and contain up to half of the stellar-mass density. However, 
hierarchical models generally predict colors for massive galaxies at z~0 that 
are too blue compared to observations (e.g., Somerville et al. 2001, Dave et al. 
2005). Some recent models suppress star formation in massive galaxies at late 
times by truncating it at some mass threshold, or by invoking strong feedback 
from AGN (e.g., Granato et al. 2001, Dave et al. 2005, Croton et al. 2006, 
Hopkins et al. 2006). The morphologies of the most optically luminous (and the 
most massive) galaxies transforms from "normal" early-type systems at z~l 
to irregular systems at z~2-3 (e.g., Papovich et al. 2005). To understand the 
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assembly of these objects we need to study the properties of massive galaxies at 
these earlier epochs. 

In these proceedings, I focus on Spitzer observations at 3-24 ^m of massive 
galaxies at z^l.5— 3 in the southern Great Observatories Origins Deep (GOODS- 
S) field. Combined with HST/ ACS and ground-based data, the Spitzer data 
provide constraints on the star-formation and AGN activity in massive galaxies 
at these epochs. I discuss the properties of the masses and IR activity in these 
galaxies, and I comment on the implications for their assembly and evolution. 

2. Stellar Masses and Star Formation in High —z Massive Galaxies 

GOODS is a multiwavelength survey of two 10'xl5' fields. The GOODS datasets 
include (along with other things) HST/ ACS and VLT/ISAAC imaging, Chandra 
X-ray observations (Giavalisco et al. 2004), and recent Spitzer imaging (see 
M. Dickinson's contribution to these proceedings). I make use of these data for 
the work described here, as well as data from Spitzer /MIPS 24 fim in this field 
from time allocated to the MIPS GTOs (e.g., Papovich et al. 2004). 

Here, I focus on so-called distant red galaxies (DRGs) selected with J—K s > 
2.3 mag (Saracco et al. 2001, Franx et al. 2003). This color criterion identifies 
galaxies at z~2-3.5 whose light is dominated by passively evolving stellar popu- 
lations older than ~250 Myr (i.e., with a strong Balmer/ 4000 A break between 
the J and -ft^-bands) , and also heavily reddened star-forming galaxies at these 
redshifts (Forster-Schreiber et al. 2004, Papovich et al. 2006). For the GOODS- 
S data, the J—K s >2.3 mag selection is approximately complete to stellar masses 
>10 n Mq for passively evolving galaxies (see Papovich et al. 2006). 

More than 50% of the DRGs have 24 /im detections with / J/ (24 / um)>50 //Jy. 
Daddi et al. (2005), Reddy et al. (2006), and Webb et al. (2006) also find similar 
24 /xm-detection rates for massive galaxies at z~2.. Interestingly, this implies 
that the majority of massive galaxies at z^2 emit strongly in the thermal IR 
- they are either actively forming stars, supermassive blackholes, or both at 
this epoch. The 24 ^m emission at z~1.5-3 probes the mid-IR (~5-10 //m), 
which broadly correlates with the total IR, Lir=L(8 — 1000/mi) (e.g., Chary & 
Elbaz 2001). Figure 1 shows the inferred Ltr for the DRGs using the Dale et al. 
(2002) models to convert the observed mid-IR to total IR luminosity. Although 
there is inherent uncertainty in this conversion, the 24 /im flux densities for the 
z~1.5-3 DRGs yield Lir~10 115 ~ 13 Lq. If attributed to star-formation, then 
this corresponds to SFRs of ^100-1000 M Q yr" 1 (e.g., Kennicutt 1998). 

Nearly all of the DRGs are detected in the deep Spitzer /IRAC data, imply- 
ing they have substantial stellar masses. In Papovich et al. (2006), we modeled 
the DRGs by comparing their ACS, ISAAC, and IRAC [3.6;um] [4.5/iin] photom- 
etry to stellar-population synthesis models (Bruzual & Chariot 2003). Although 
the modeling loosely constrains the ages, dust content, and star-formation his- 
tories of the DRGs, it provides relatively robust estimates of the galaxies' stellar 
masses, which are plotted in figure 1 as a function of redshift. Typical uncer- 
tainties for the stellar masses for the full DRG sample are 0.1-0.3 dex. 

Figure 2 shows the specific SFRs ($f/M) derived from the masses and 
SFRs for the DRGs, where the SFRs are derived from the summed UV and IR 
emission. The figure also shows the specific SFRs for lower redshift galaxies 
from COMBO-17 (Wolf et al. 2003), which overlaps with the GTO 24 /im data. 



Spitzer Observations of Red Galaxies at High Redshifts 



10 1 



O DRG f„(24^m) ^ 50/xJy 
k DRG X-ray detection 
□ HDF-N 24-fj.m sample 




10 12 






X O 


° 




h 






°.§ 

*9 □ 





• X 
• X 


y 







- -cm - • -ffl -« or 


&-*"• 






10" 




•\°.° 






- 






O ^ 


"V*' «< 


a 










o, v^' 


□ ° ° 








• 




if" 


□ 






10'° 


/ 




• 






1: 










■ DRG, no 


24 fj,m 










o DRG, 24 


fj,m 












* X-ray DRG 
















10" 












i 



Figure 1. LEFT: Total IR luminosities, Lir=L(8 — 1000), of galaxies in- 
ferred from their observed 24 /im emission. Symbol definitions are inset. The 
solid line denotes the 24 /zm 50% completeness limit of the GTO data. The 
estimated systematic error is «0.5 dex.. RIGHT: Stellar masses of galaxies 
inferred by fitting models to the galaxies' rest-frame UV-to-near-IR data. 
The inset bars show the mean errors as a function of mass. The short-dashed 
line shows the characteristic present-day stellar mass (Cole et al. 2001); the 
long-dashed line shows the stellar mass limit for a passively evolving stellar 
population formed at z^oo with K s =23.2 mag. 

The massive galaxies with _M>10 n A4© at 1.5<z<3 have high specific SFRs, 
\I//.A/f~0.2-10 Gyr -1 (excluding X-ray sources). In contrast, at z<^ 0.75 galaxies 
with _M>10 n M & have much lower specific SFRs, ^/A4~0.1-l Gyr -1 . 

In Papovich et al. (2006), we defined the integrated specific SFR as the 
ratio of the sum of the SFRs, ^i, to the sum of their stellar masses, M.- L , 
T=J2i ^?i/J2i Mi, summed over all i galaxies. Figure 3 shows the integrated 
specific SFRs for DRGs at s~l. 5-3.0 and COMBO-17 at z ~ 0.4 and 0.7, all 
with A4 > 10 11 A4q. The error box indicates the affect of assumptions in the 
SFRs and AGN activity in the DRGs (see further discussion below, and in Pa- 
povich et al. 2006). The integrated specific SFR in galaxies with A / (>10 11 Mq 
declines by more than an order of magnitude from z~1.5-3 to z<^0.7. 

The curves in figure 2 show the specific SFR integrated over all galaxies (not 
just the most massive); this is the ratio of the cosmic SFR density to its integral, 
p*/ J p* dt. Although there is a decrease in the global specific SFR with decreas- 
ing redshift, the evolution in the integrated specific SFR in massive galaxies 
is accelerated. The implication is that at z^>1.5, massive galaxies are rapidly 
forming their stars, whereas by z<^1.5 the specific SFRs of massive galaxies drops 
rapidly, and lower-mass galaxies dominate the cosmic SFR density . 

3. The Contribution of AGN to the mid IR Emission 



Many (~15%) of the massive galaxies at z~1.5-3 are detected in the deep X-ray 
data, and these objects tend to have high inferred IR luminosities and specific 
SFRs (see figs. 1 and 2). Indeed, the fraction of X-ray detected DRGs rises 
with increasing luminosity (>50% at Ltr^IO 13 L©), similar to the detection 
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Figure 2. Specific SFRs as a function of galaxy stellar mass. The Left 
panel shows the DRG and HDF-N galaxies with 1.5 <z< 3, with symbols 
the same as in figure 1. The Right panel shows lower-redshift galaxies from 
COMBO-17. Open circles show COMBO-17 galaxies with 24 /xm detections, 
small filled symbols show upper limits for galaxies undetected at 24 /jm. 



rate of sub-mm galaxies (Alexander et al. 2005). At these redshifts the X-ray 
fluxes imply the presence of an AGN with Lx^lO 42 erg s _1 . Many authors 
are also finding AGN candidates based on rest-frame near-IR colors, where 
the AGN is presumably heavily obscured by gas and dust so that it is missed 
in deep (^1 Msec) X-ray surveys (e.g., Donley et al. 2005, Stern et al. 2005, 
Alonso-Herrero et al. 2006, Barmby et al. 2006). Approximately 10% of the 
X-ray undetected DRGs have ACS-through-IRAC colors consistent with dust- 
enshrouded AGN. Combined with the 15% of DRGs detected in the X-rays, up 
to 25% of the DRG population host AGN (see also Papovich et al. 2006). 

If AGN contribute to the observed 24 /im emission in galaxies at z~1.5-3, 
then they can affect the inferred IR luminosities. For example, using an IR 
template for Mrk 231 instead of a star-forming galaxy with Ltr^IO 13 L would 
reduce the IR luminosity for z~1.5-3 galaxies by a factor of ~2-5. To limit this 
effect on the evolution of the integrated specific SFRs, the error box in figure 3 
shows how the result changes if for galaxies with putative AGN we set the SFR 
to zero. The high AGN occurrence in DRGs suggests that massive galaxies at 
z~1.5-3 simultaneously form stars and grow supermassive black holes. 

However, even when AGN are present it is unclear whether star-formation 
or the AGN activity domninates the bolometric IR luminosity. Although ~80% 
of sub-mm galaxies have X-ray detections (Alexander et al. 2005), the IR to 
X-ray luminosity ratios are up to an order of magnitude higher than what is 
expected for AGN alone. This suggests that both star-formation and AGN 
contribute to the bolometric emission. Similarly, in the near-IR spectrum of 
a z~2.5 DRG, Frayer et al. (2003) find that the galaxy nucleus has a low Ha 
to [Nil] flux ratio consistent with ionization from an AGN. However, Ha is 
spatially resolved in their spectrum and the [Nil] line strength drops off in the 
off-nucleus spectrum. This implies extended star-formation beyond the nucleus, 
which presumably contributes to the inferred IR luminosity. Both AGN and 
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Figure 3. Evolution of the integrated specific SFR, the ratio of the total 
SFR to the total stellar mass (from Papovich et al. 2006). The curves show the 
expected evolution from the global SFR density (solid lines, Cole et al. 2001, 
thick line includes correction for dust extinction; dashed line, Hernquist & 
Springel 2003). Data points show results for galaxies with >10 n M-q- Filled 
circle corresponds to the DRGs; filled diamonds correspond to the COMBO- 
17 galaxies. The inset error bar shows an estimate on the systematics. 

star-formation occur simultaneously in high-redshift IR-detected galaxies and 
both probably contribute to the IR emission. 



4. Summary 

In this contribution, I discussed star-formation and AGN activity in massive 
galaxies (;>10 n M.®) at z~l-3 using observations from Spitzer at 3-24 /xm. The 
majority (^50%) of these objects have /„(24/im)>50 /iJy, which if attributed 
to star formation implies SFRs of ;>100 Mq yr _1 . Galaxies at z~1.5-3 with 
.A/f >10 n M.q have specific SFRs equal to or exceeding the global average value. 
In contrast, massive galaxies at z~0.3-0.75 have specific SFRs less than the 
global average, and more than lOx lower than that at z~1.5-3. By z-^l.b 
massive galaxies have formed most of their stellar mass, and lower-mass galaxies 
dominate the SFR density. At the same time, as many as 25% of the massive 
galaxies at z^l.5 host AGN. The high AGN occurrence at z~ 1.5-3 provides 
evidence that massive galaxies are simultaneously forming stars and growing 
supermassive black holes. This may provide the impetus for the present-day 
black-hole-bulge-mass relation and/or provide the feedback necessary to squelch 
star-formation in such galaxies, moving them onto the red sequence. 
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Lastly, while on average high-redshift galaxies have IR spectral energy dis- 
tributions consistent with local templates, individually there remains signifi- 
cant uncertainty. Future work is needed to understand the distribution between 
the mid-IR (rest-frame 5-15 /jm) emission and total IR luminosity in z~1.5-3 
galaxies. Herschel Space Observatory (and eventually SAFIR) will mitigate this 
problem by measuring the far-IR emission of distant galaxies directly. 
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